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detonated  above  or  below  VVT.  Results  from  the  analysis  of  17  N  TS  shots  recorded  at  the  RSTN  sta¬ 
tion,  RSSD  are  generally  not  as  well  defined  but  appear  similar  to  those  from  the  JAS  recordings.  The 
largest  difference  between  die  below  and  above  WT  shots  of  a  given  yield  lies  in  their  low-  frequency 
spectral  level  and,  to  a  slightly  smaller  extern,  in  their  rnb  values. 

Examination  of  mh  versus  yield  from  large  data  bases  of  NTS  shots  recorded  at  regional  and 
teleseismic  distances  also  shows  that,  for  a  fixed  yield,  shots  below  WT  have  significantly  larger  mb 
than  those  detonated  above  WT.  Gas-filled  porosity  of  die  shot  medium  appears  to  be  die  dominant 
reason  for  the  large  magnitude  bias  between  below  and  above  WT  shots,  l'or  shots  detonated  near  die 
water  table,  explosion  cavity  and  its  proximity  to  WT  can  also  significantly  inlluence  the  source  cou¬ 
pling.  Below  and  above  WT  shots  may  be  distinguished  by  comparing  their  Rn  and  Lg  phases  at  com¬ 
mon  recording  stations. 
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S  LAI  MARY 


A  comparison  of  the  spectra  of  nuclear  shots  detonated  below  and  above  the  water  table  (WT)  ts 
made  by  analysis  of  mostly  the  regional  phase  Pn  recorded  at  two  common  stations.  Parameters 
derived  from  each  record  include  estimates  of  corner  frequency  f  and  the  low -frequence  spec'ral  level 
A0.  after  the  spectra  have  been  corrected  tor  the  effects  of  spatial  attenuation.  Tne  tine  domain  meas¬ 
urements  of  A/T.  a  measure  of  mh,  are  also  made  on  each  observed  record  Results  from  22  Nevada 
Test  Site  (NTS)  shots  recorded  at  the  DWWSSN  digital  station,  JAS.  indicate  that  plots  of/,  vet sus  A„ 
show  some  separation  between  the  below  and  above  VVT  populations  whereas  plots  of  /  versus  A  T 
show  no  significant  separation.  I  or  a  fixed  yield,  f.  for  shots  below  and  above  WT  a.s  ‘ariv  indis¬ 
tinguishable  whereas  both  A0  and  A/I'  are  considerably  larger  for  shots  below  WT  than  for  those  above 
WT,  presumably  due  »•>  improved  source  coupling  of  shots  m  water  saturated  environment.  A  com 
partson  of  yields  estimator  from  A0,  AT,  A0  combined  with  f .  and  AT  combined  w  ith  /,  demonstrates 
that  substantially  more  precise  yield  determinations  can  be  made  h  o  Known  whether  the  shot  was 
detonated  above  or  below  WT.  Results  from  the  analysis  ot  17  NTS  shots  recorded  at  rite  RSTN  sta¬ 
tion.  RSS1)  are  generally  not  as  well  defined  but  apjx'.ir  similai  to  those  iiom  the  JAS  recordings.  The 
largest  difference  between  the  bel  >w  and  above  WT  shots  of  a  given  v  icld  lies  in  thou  low  -  frequenev 
spectra!  level  and,  to  a  siightly  sma.ler  extent,  in  theii  rnh  values 

Examination  ot  mh  versus  yield  ironi  large  data  bases  ot  N1S  shots  recorderl  at  regional  and 
teleseismic  distances  also  shows  that,  for  a  fixed  yield,  shots  below  WT  have  signilicantk  larger  m, 
than  those  detonated  above  WT  (ias  Idled  porosity  of  the  shot  medium  appeals  to  be  the  dommant 
reason  tor  the  latge  magnitude  bias  between  tvlnw  and  above  W!  shorn  I'oi  slots  detonated  neat  the 
water  table,  explosion  cav  itv  and  its  pioximitv  to  Wl  can  al-o  s;  :n: ncautl v  mtluence  tile  soutce  con- 
pi 1 1 1 Below  ap  i  shove  WT  shots  mav  tv  distinguished  b\  rompaiing  then  Pn  and  1  .:  phases  at  iom 
m>  m  lev  i udme  stale  >:is 
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INTRODUCTION 


The  geological  environment  in  which  a  nuclear  explosive  is  emplaced  has  a  major 
influence  on  the  strength  of  seismic  signals.  For  example.  Werth  and  Herbst  (196.')  tound  that 
for  a  fixed  yield,  the  first  half-cycle  amplitudes  for  granite  and  salt  are  nearly  equal,  for  tuff 
the  amplitude  is  about  one-half  the  amplitudes  for  granite  and  salt,  and  for  dry  alluvium  the 
amplitude  is  about  one-tenth  the  amplitudes  for  granite  and  salt.  These  trends  were  confirmed 
by  Springer  (1966)  who  further  suggested  that  dry  porosity  is  the  parameter  responsible  for 
such  wide  variations  in  seismic  coupling.  Murphy  (1981)  reviewed  available  observed 
mb/yield  and  free-field  data  from  contained  explosions  in  a  variety  of  source  media  and  con¬ 
cluded  that  the  coupling  into  the  teleseismic  P  wave  transmission  path  is  about  the  same  for 
shots  in  granite,  shale,  and  wet  tuff-rhyolite  emplacement  media,  whereas  explosions  in  dry 
unconsolidated  material  couple  significantly  less  efficiently 

Worth  and  Herbst  (1963)  noted  that  explosions  in  higher  velocity  rock  arc  significantly 
richer  in  higher  frequencies.  Yet  their  data  for  5  kt  shots  in  alluvium  and  tuff  (Figure  3, 
Werth  and  Herbst,  1963)  indicate  more  high  frequencies  for  alluvium  than  for  tuff  (or  higher 
corner  frequency  fc  for  alluvium  than  for  tuff),  together  with  amplitudes  in  tuff  about  5  times 
larger  than  in  alluvium.  Non-linear  finite  difference  calculations  (Ruche.  1982)  indicate  that 
the  most  important  characteristic  of  NTS  tuff  is  its  porosity,  which  can  account  tor  the  large 
observed  diffetenccx  in  the  seismic  coupling  of  events  above  and  below  the  water  table.  A 
comparison  of  source  functions  tor  Yucca  Flat  wet  and  dry  tut!  (Rache.  1982;  Figure  lit 
shows  a  larger  low-frequency  amplitude  and  lower  tv  for  wet  tuff  in  comparison  to  dip,  tuff, 
when  vield  values  are  the  same  Fvemden  <  t  al.  il9X0)  also  show  similar  strong  dependence 
of  explosion  I’  wave  spcUra  on  the  mateiial  poverties  of  the  coupling  media  (see  especially 


I 


their  Figure  8).  Gupta  and  Blandford  (1987)  examined  teicseismic  I’  hum  shots  in  saturated 
(below  water  table)  media  and  found  spectral  dillerenees  to  depend  on  the  shot  medium  vein 
city.  It  seems  therefore  that  both  amplitudes  and  spectral  shapes  are  significantly  influenced 
by  the  geological  environment  of  the  shot  point. 

In  this  study,  the  amplitude  and  sjrectra  of  the  regional  phase  Pn  are  compared  for  shuts 
detonated  below  and  above  the  water  table.  The  three  parameters  estimated  from  each  record 
for  this  purpose  are  comer  frequency  f.,  low-frequency  spectral  level  A0,  and  the  nme-domam 
measurement  AH'  (a  measure  of  nij,).  Short  period,  vertical -component  data  from  two  sta¬ 
tions,  JAS  (with  records  from  22  shots)  and  RSSI)  (with  data  from  17  shots)  are  analyzed 
I  he  results  indicate  that,  for  a  gi\ en  yield,  botli  A,,  and  mb  arc  significantly  larger  for  shots 
below  \\  1  than  tor  those  above  W  I ,  but  the  f c  values  are  not  distinguishable.  These  results 
are  in  agreement  with  those  derived  from  substantially  larger  data  sets  of  both  regional  and 
t  e  1  e  s  e  i  s  i  n  i  e  measurement" 


ANALYSIS  OF  REGIONAL  DATA  RECORDED  AT  JAMESTOWN  tJAS) 


Spectral  and  Time  domain  Measurements: 

We  analyzed  digital  data  from  22  NTS  shots  well  recorded  at  20  samples/sec  at  the 
DWWSSN  station  Jamestown,  JAS  at  epicentral  distances  of  about  350  km.  Table  1  lists 
these  22  explosions  along  with  pertinent  information  such  as  the  shot  medium,  magnitude  mb, 
and  whether  the  shot  point  is  above  or  below  the  water  table  Out  o!  22  shots,  onlv  9  are 
below  WT;  19  shots  are  from  the  Yucca  Hats  and  only  3  (SHRI’A.  HARZHR.  and  SAl.l'T) 
are  from  the  Pahute  Mesa  regions  of  the  NTS 


TABU-  1 


NTS  SHOTS  RECORDED  AT  JAS 


Shot 

G  as 

Shoi 

I  >epth 

Nj. 

Date 

Name 

Medium 

mh 

Porositv 

Depth 

of  WT 

Symbol* 

(vol  %) 

(m) 

( ml 

1 

14  Nov  1980 

DAUPHIN 

hjfT 

4.1 

17.0 

320 

580 

YO 

17  Dec  1980 

SERPA 

tuff 

5  1 

10  0 

573 

627 

PO 

3 

29  May  1981 

aligou: 

tuft 

4  2 

8.0 

320 

605 

YO 

4 

06  Jun  1981 

HAR7.F.R 

tuff 

5.5 

3.0 

637 

668 

PO 

5 

27  Aug  1981 

1S1.AY 

tuff 

- 

23.0 

294 

567 

YO 

6 

01  Oct  1981 

PAL1ZA 

tuff 

4.9 

6.0 

472 

570 

YO 

7 

1!  Nov  1981 

T1LC1 

alluvium 

4  8 

10.1 

445 

464 

YO 

8 

12  Nov  1981 

ROUSANNE 

tuff 

5  3 

•2  0 

518 

495 

YX 

9 

03  Dec  1981 

AKA  VI 

tuff 

4.6 

14.3 

494 

580 

YO 

10 

28  Jan  1982 

JORNADA 

tuff 

5.9 

0.0 

640 

507 

YX 

1 1 

29  Sep  1982 

BORREGO 

tuff 

00 

56-1 

501 

YX 

i: 

26  May  1983 

FAHADA 

tuff 

4.4 

12  0 

384 

600 

YO 

13 

09  Jun  1983 

DANABLIJ 

alluvium 

4.5 

12.5 

320 

584 

YO 

14 

22  Sep  1987 

TECH  AIK) 

tuft 

0  0 

573 

500 

YX 

15 

31  May  1984 

CAPROCK 

tuff 

5.8 

0.0 

(>00 

500 

YX 

lo 

20  Jun  1984 

Dl'ORO 

tuff 

4.6 

14.0 

381 

480 

YO 

r 

02  Aug  1984 

CORRLO 

tuff 

4.7 

1  3.0 

375 

470 

YO 

18 

15  Mar  1985 

VAUGHN 

tuff 

4,8 

110 

427 

498 

YO 

19 

12  Jun  1985 

SAIT  'I 

thyolite 

5  5 

4.0 

698 

622 

PX 

.10 

05  Dec  1985 

KINTBOTO 

tuff 

5.7 

0.0 

579 

488 

YX 

:i 

22  Mar  1986 

GLENCOE 

tuff 

5.1 

0  0 

f>]0 

822 

YX 

“V  ■> 

14  Nov  1986 

GASCON 

tuff 

5  8 

0  0 

594 

51  IS 

YX 
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Mesa  and  3  ucc, 

i  Mats  regions 
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X  below  (he  water  table  shots 
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et location  Experiment 
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ot.  KEARSAKCiE  (17  Aug 
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7  29" 
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m  Figure  1 

1  he  JVE 

shot  was  tec 

orded  at  the 

new  DWWNSN  station,  CM!!,  ahi 

11;  1(1  kr 
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available  data 
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ow  33  r, 

has  mH 

,5 . 9  w  h  1  e  h 

suggests 

a  yield  value 

close  to  1st) 

k  t  ( sec .  c  r 

Mac! 

he,  1 9S2 ) 

I  lie  .17 

.  E  shut. 

known,  to  be 

above  33 

1  .  also  had  a 

\  icid  ot  pioh. 

ibis  elosc  'a ' 

1  so  kt  althoug 

h  its  in,, 

is  onlv  .8  4  The  huge 

ditlerence  m  the  1’n  amplitudes  loi  the  two  shots  is  most  probably  due  to  the  large  ditlerences 


J 


JORNADA  (JAS)-BWT 


and  KEARSARGE  (above  water  table)  at  CMB  (located  within  10  km  of  JAS)  v  uii  2  sec  of 
noise  before  the  onset  of  Pn,  with  two  different  time  scales.  Note  the  large  difference  in  the 
Pn  amplitudes.  Pg  is  clipped  for  both  events  but  Pn/Pg  amplitude  ratio  is  mud  mealier  for 
KEARSARGE  than  for  JORNADA. 


in  coupling  expected  for  shots  below  and  above  the  WT.  It  is  interesting  to  note  that  although 
Pg  is  clipped  for  both  shots,  the  Pn/Pg  amplitude  ratio  is  much  smaller  for  KEARSARGE 
than  for  JORNADA,  in  fact,  a  comparison  of  the  later,  unclipped  portion  of  Pg  (or  Pg  "coda") 
amplitudes  appears  to  suggest  that  Pg  may  not  show  the  effect  of  the  water  table. 

A  window  length  of  6.4  sec  (’28  points)  with  10%  cosine  taper  is  used  for  spectral 
analysis  of  Pn.  The  spectra  are  corrected  for  noise  by  computing  spectra  of  a  window  of  noise 
prior  to  the  onset  of  Pn  and  subtracting  the  power  in  the  noise  from  the  power  in  the  observed 
signal  for  each  discrete  frequency.  Examples  of  such  spectra  of  the  observed  signal  and 
noise,  up  to  the  Nyquist  frequency  of  10  Hz,  are  shown  in  Figure  2  for  the  same  'wo  records 
as  in  Figure  1.  The  spectra  are  without  smoothing  and  not  corrected  for  instrumental  response. 
The  signal-io-noise  (S/N)  ratio  is  good  for  almost  the  entire  frequency  range.  The  spectrai 
ratio  IvEARSARGE/JORNADA,  corrected  for  noise,  and  shown  in  Figure  3,  has  a  mean  slope 
of  0.1 12  per  Hz  (with  associated  standard  deviation,  SD  of  0.004)  over  the  frequency  range  of 
0.2  to  6.0  Hz.  Results  for  the  frequency  range  of  0.2  to  10.0  Hz  are  also  shown  in  the  figure 
The  mean  values  of  the  spectral  ratio  (in  magnitude  units),  integrated  over  the  two  frequency 
ranges  by  using  the  mean  values  of  slope  and  intercept,  denoted  by  AV  1NTEG  in  Figure  3, 
indicate  the  average  spectral  amplitudes  in  KEARSARGE  to  be  nearly  half  that  for  JOR¬ 
NADA. 

The  far-held  source  spectra  may  generally  be  characterized  by  three  independent  parame¬ 
ters:  the  low-frequency  spectral  level  A0,  proportional  to  seismic  moment;  the  comer  fre¬ 
quency  fc,  and  the  power  of  the  high-frequency  asymptote  (Aki  and  Richards,  1980)  The  last 
parameter  or  the  value  of  n  in  f  n  may  be  assumed  to  be  2  for  explosions  (Evemden  ct  ol  , 
1986).  For  explosions  in  same  or  similar  media,  A0  is  a  good  measure  of  the  explosion  yield 
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Figure  2.  Vertical-component  displacement  amplitude  spectra  (symbol  +)  of  Pn,  uncorrected 
for  instrument  response  and  based  on  6.4  sec  long  window  with  10%  cosine  taper,  on  JAS 
and  CMB  records  of  JORNADA  and  KFARSARGE,  respectively.  Spectra  of  an  equal  win¬ 
dow  length  of  noise  (symb^1  N  ure  also  included.  S/N  appears  good  up  to  the  Nyquist  fre¬ 
quency  of  10  Hz. 


KEARSARGE  (CMB)/  JORNADA  (JAS) 
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Oo  on  ~~2.  00  00  60  00  ii.  00 

FREQUENCY  (Hz) 


FREQ  BAND  SLOPE  S.D.  AV  INTEG 
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Figure  3.  Spectral  ratio  KEARSARGE/JORNADA  of  Pn,  corrected  for  noise.  Points  for  which 
S/N  power  ratio  is  less  than  1.5  are  not  plotted.  The  dashed  line  shows  the  mean  least  squares 
slope  over  the  frequency  range  of  0.2  to  6.0  Hz.  Mean  slope,  with  associated  standard  devia¬ 
tion,  SD  and  mean  value  of  log  amplitude  ratio  (AV  INTEG)  are  indicated  for  two  frequency 
ranges. 
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(Evemden  etal.,  1986;  McLaughlin  etal.,  1987).  A  spectral  measure  of  A0  is  made  by 
correcting  for  noise,  removing  the  instrument  response,  applying  the  appropriate  spatial 
attenuation  correction  and  then  least-squares  fitting  the  log-amplitude  spectra  to  a  model  spec¬ 
trum  with  a  fall-off  rate  of  f  "2.  Only  those  data  points  that  have  S/N  larger  than  a  specified 
value  are  used  in  the  least-squares  regression.  The  amplitude  spectrum  is  assumed  to  have  the 
form  (Brune,  1970) 


A(f)  = 


1  +  (f/fc)2 


Bakun  and  Johnson  (1970)  analyzed  earthquake  and  NTS  explosion  data  recorded  at  JAS 
and  found  Q  to  be  about  400  for  Pn  at  regional  distances.  This  corresponds  to  t*  of  about  0.1 
for  NTS  shots  recorded  at  JAS.  Applying  this  attenuation  correction,  spectral  fits  for  JOR¬ 
NADA  and  KEARSARGE,  based  on  S/N  power  threshold  of  1.5,  are  shown  in  Figure  4 
which  clearly  demonstrates  the  latter  shot  to  have  significantly  lower  A0  and  higher  fc  than 
those  for  the  former  shot.  The  frequency  range  used  for  spectral  fits  was  0.2  to  6.0  Hz  since 


most  shots  showed  good  S/N  only  in  this  frequency  band. 


A  similar  procedure  was  applied  to  all  22  shots  to  obtain  values  of  A0  and  fc,  which  are 
plotted  on  log-log  scale  in  Figure  5a  in  which  the  explosion  numbers  correspond  to  those  in 
Table  I.  Below  WT  (denoted  by  X)  and  above  WT  (marked  by  circles)  shot  populations  do 
not  appear  to  show  an  obvious  separation.  The  mean  slope  and  intercept  values  are  -0.245 
(0.051)  and  0.654  (0.104),  respectively,  where  the  quantities  within  parentheses  denote  one 
SD  value.  Data  from  the  9  below  WT  shots  and  the  13  above  WT  shots  are  plotted  in  Figures 
5b  and  5c,  respectively.  Least  squares  linear  regressions  of  data  in  these  two  figures  yielded 
mean  slopes  of  -0.302  (0.061)  and  -0.283  (0.104),  respectively.  Both  of  these  slope  values  are 
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results  are  listed. 


[V  i 


Figure  5a.  fc  versus  A0  for  9  below  water  table  (denoted  by  X)  and  13  above  water  table 
(marked  by  circles)  shots  numbered  as  in  Table  1  and  based  on  t*  =0.1.  The  least  squares 
linear  regression  through  all  22  data  points  (continuous  line)  has  correlation  coefficient  of 
0.731,  slope  of  -0.245,  and  intercept  of  0.654.  One  SD  variation  in  the  intercept  value  is 
shown  by  the  dashed  lines. 
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Figure  5b.  Similar  to  Figure  5a  for  only  9  below  water  table  shots.  Forcing  a  slope  of  -1/3, 
the  least  squares  linear  regression  (continuous  line)  has  an  intercept  of  0.926. 
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Figure  5c.  Similar  to  Figure  5a  for  only  13  above  water  table  shots.  Forcing  a  slope  of  1/3, 
the  least  squares  linear  regression  (continuous  line)  has  an  intercept  of  0.754. 
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within  one  SD  of  the  theoretically  expected  value  of  -1/3.  We  therefore  forced  a  slope  of  -1/3 
through  the  data  points  in  Figures  5b  and  5c  and  the  least  squares  regressions  gave  mean 
intercept  values  of  0.926  (0.044)  and  0.754  (0.045),  respectively.  Thus  the  below  and  above 
WT  populations  are  separated  by  0.172  (0.063).  In  other  words,  for  the  same  low-frequency 
spectral  level,  shots  below  WT  have,  on  the  average,  larger  (about  a  factor  of  1.5)  fc  than 
those  above  WT. 

Determination  of  fc  is  obviously  uependent  on  the  attenuation  relationship  assumed  in 
correcting  the  observed  amplitudes.  A  different  attenuation  correction  is  therefore  tried  in 
order  to  examine  how  the  results  pertaining  to  below  and  above  WT  shot  populations  are 
affected.  Analyzing  Pn  spectra  of  a  large  number  of  NTS  shots  recorded  at  regional  dis¬ 
tances,  Taylor  et  al.  (1987)  suggested  the  following  frequency-dependent  attenuation  relation¬ 
ship: 

Q  =  250  f  0  6  (2) 

This  correction  is  applied  to  the  data  by  assuming  a  propagation  velocity  of  8  km/sec  for  Pn. 
The  resulting  plots  of  A0  versus  fc  are  shown  in  Figures  6a,  6b,  and  6c  for  all  (22)  shots,  9 
shots  below  WT,  and  13  shots  above  WT,  respectively.  The  mean  slope  and  intercept  values 
for  all  shots  (Figure  6a)  are  -0.237  (0.042)  and  0.497  (0.098),  respectively.  The  mean  slope 
values  for  the  data  in  Figures  6b  and  6c  are  -0.276  (0.055)  and  -0.272  (0.087),  respectively; 
both  of  these  slope  values  are  within  one  SD  of  the  theoretically  expected  value  of  -1/3.  A 
slope  of  -1/3  is  therefore  forced  through  the  data  points  in  Figures  6b  and  6c  and  the  least 
squares  regressions  give  mean  intercept  values  of  0.802  (0.044)  and  0.648  (0.039),  respec¬ 
tively.  Thus  the  below  and  above  WT  populations  are  separated  by  0  154  (0.059).  This  result 
is  very  neaiiy  the  same  as  obtained  earlier  based  on  t*  =  0. 1.  It  means  that  the  separation 
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Figure  6a.  Similar  to  Figure  5a  but  based  on  frequency-dependent  Q.  The  least  squares  linear 
regression  (continuous  line)  has  correlation  coefficient  of  0.782,  slope  of  -0.237,  and  intercept 
of  0.497. 
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between  below  and  above  WT  populations  is  not  much  influenced  by  the  choice  ot  attenuat'.or. 
model  and  is  therefore  significant. 

Note  that  as  far  as  the  distinction  between  beiow  and  above  WT  shots  is  concerned,  the 
results  based  on  t*  =0.1  (Figures  5a,  b,  c)  and  those  based  on  frequency-dependent  Q  (Fig¬ 
ures  6a,  b,  c)  are  nearly  identical;  the  relative  distribution  of  various  data  points  is  almost  the 
same  for  the  two  sets  of  figures  and  only  the  mean  intercept  values  are  different.  Also  the 
separation  between  the  below  and  above  WT  populations  is  nearly  independent  of  the  attenua¬ 
tion  model  used  for  correcting  Pn  amplitudes  For  this  reason,  subsequent  analyses  of  data 
are  restricted  to  the  use  of  only  one  attenuation  model,  viz.  t*  =0.1. 

Time  domain  measurements  of  AT  are  also  made  on  all  22  JAS  records.  Plots  of  fc 
(derived  by  assuming  t*  0  1)  versus  log  AT  aie  shown  in  Figures  7a,  7b,  and  7c  for  all 
(22)  shots,  below  WT  (9)  shots,  and  above  WT  (15)  shots,  respectively.  The  mean  siupc 
values  for  the  data  in  these  three  figures  are  0  220  (0080),  -0  296  (0  110),  and  -0.149 
(0.164).  respectively  When  the  mean  slope  derived  from  the  entire  data  set  (viz.  -0.220)  is 
forced  through  the  data  points  in  f  igures  7b  and  7c.  the  mean  intercept  values  are  0.714 
(0  067)  and  0.654  (0.056),  respectively  Thus,  for  a  given  AT  or  mh  value,  f.  values  for  shots 
beiow  and  above  WT  are  not  significantly  different 

Plots  ot  t  versus  vield  (on  log  log  scale)  be  all  shots  showed  mean  slope  and  intercept 
values  of  o  249  (0070)  and  0  541  (0  109),  respectively  For  below  WT  and  above  WT  shots, 
t’ne  mean  slopes  were  ()  24’7  (OO'T)  and  0  24')  (0  118),  respectively.  Both  of  these  slo|v 
values  aie  within  one  SI)  of  the  theoretically  expected  value  of  15  A  slope  of  -1/5  was 
fherctorc  forced  through  the  data  points  and  the  least  squares  regressions  gave  mean  intercept 


IX 


log  A/T 

bigure  7a.  fc  versus  log  ATT  for  9  below  water  table  (denoted  by  X)  and  13  abose  water  ta.ble 
(marked  by  circles)  shots  based  on  t*  =  0.1.  The  least  squares  linear  regression  through  all  32 
data  points  (continuous  line)  has  correlation  coefficient  of  0.524,  slope  of  -0  320,  and  intercept 
of  0.678.  One  SD  variation  in  the  intercept  value  is  shown  by  the  clashed  lines 
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values  of  0.663  (0.069)  and  0.662  (0.051),  respectively.  Thus  there  is  no  separation  between 
the  below  and  above  WT  populations.  In  other  words,  for  a  fixed  yield,  fc  values  for  shots 
below  and  above  WT  are  not  significantly  different. 

Plots  of  A0  versus  yield  (on  log-log  scale)  for  all  (22)  shots,  below  WT  (9)  shots,  and 
above  WT  (13)  shots  gave  mean  slopes  of  1.061  (0.124),  0.969  (0.118),  and  1.094  (0.119), 
respectively.  All  of  these  slope  values  are  within  one  SD  of  the  theoretically  expected  value 
of  1.  When  a  slope  of  1  was  forced  through  the  data  points,  the  least  squares  regressions  give 
mean  intercept  values  of  0.789  (0.083)  and  0.276  (0.045),  respectively.  Thus  the  below  and 
above  WT  populations  are  separated  by  0.513  (0.095).  In  other  words,  for  the  same  value  of 
yield,  A0  for  a  shot  below  WT  is,  on  the  average,  significantly  larger  (about  a  factor  of  3) 
than  that  for  a  shot  above  WT. 

Plots  of  log  AT  versus  yield  for  all  shots,  below  WT  shots,  and  above  WT  shots  gave 
mean  slopes  of  0.875  (0.086),  0.829  (0.048),  and  0.846  (0.1 1 1),  respectively.  When  the  mean 
slope  derived  from  the  entire  data  set  (viz.  0.875)  was  forced  through  the  data  points,  the 
mean  intercept  values  were  1.206  (0.036)  and  0.844  (0.041),  respectively.  Thus,  the  below  and 
above  WT  populations  are  separated  by  0.362  (0.0542).  In  other  words,  for  a  fixed  yield,  log 
AT  or  mb  for  a  shot  below  WT  is  significantly  larger  (about  0.36  magnitude  unit)  than  that 
for  a  shot  above  WT. 

Yield  Estimation: 

In  order  to  assess  the  effect  of  the  water  table  on  yield  estimates,  we  computed  estimated 
yields  (Y)  based  on  A0,  AT,  A0  and  ft,  and  AT  and  T  The  coefficients  used  to  compute  the 
estimates  are  first  determined  through  linear  or  multilinear  regressions,  with  the  known  yield 
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(Y)  as  the  dependent  variable: 


log  Y  =  a,  +  b,  log  A0 

(3a) 

log  Y  =  a2  +  b2  log  A/T 

(3b) 

log  Y  =  a3  +  b3  log  A0  +  c3  log  fc 

(3c) 

log  Y  =  a4  +  b4  log  A/T  +  c4  log  fc 

(3d) 

Regression  coefficients  as  well  as  the  estimated  yields  (Y)  for  each  shot  (along  with  its  associ¬ 
ated  SD)  are  computed  for  9  events  below  the  water  table  (abwt,  b,bwl,  cbwl),  for  13  events 
above  the  water  table  (a*wt,  b*wt,  ciawl),  and  for  all  (22  shots)  of  the  data  together 
(a,al1,  b3'1,  c311).  The  regression  results  based  on  log  Y  as  the  dependent  variable  are  listed  in 
Table  2  in  which  the  SD  of  the  estimated  yield  (in  magnitude  units)  is  denoted  by  s. 

TABLE  2 


LINEAR  REGRESSION  COEFFICIENTS  FOR  JAS  DATA 


Data 

Set 

d* 

Constant 

a 

Independent 

Variable 

Coefficient 

b 

Independent 

Variable 

Coefficient 

c 

s 

ALL(22) 

20 

0.014 

0.740 

0.276 

BWT(9) 

7 

-0.638 

log  A0 

0.935 

0.261 

AWT(13) 

1 1 

0.040 

0.809 

0.141 

ALL(22) 

20 

-0.718 

0.958 

0.239 

BWT(9) 

7 

-1.469 

log  A/T 

1.178 

0.129 

AWT(13) 

1 1 

-0.663 

0.995 

0.165 

ALL(22) 

19 

-0.081 

0.776 

0.145 

0.282 

BWT(9) 

6 

-2.095 

log  A0 

1.451 

log  fc 

1.707 

0.115 

AWT(13) 

10 

-0.211 

0.915 

0.374 

0.132 

ALL(22) 

19 

-0  391 

0.852 

-0.482 

0.224 

BWT(9) 

6 

-1.446 

log  A/T 

1.171 

log  fc 

-0.025 

0.139 

A  WT(  1 3) 

10 

-0.414 

0.923 

-0.490 

0.136 

*  d  denotes  degrees  of  freedom 
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In  order  to  assess  the  improvement  in  yield  estimates  obtained  by  considering  events 
above  the  water  table  separately  from  events  below  the  water  table,  we  computed  two  sets  of 
estimated  yields  using  the  coefficients  in  Table  2.  The  first  set  is  computed  using  the  regres¬ 
sion  coefficients  obtained  for  all  of  the  data  together,  e.g., 

log  Y j  =  af11  +  b,aIi  Gog  A0)j  (4) 

The  second  set  of  yield  estimates  is  computed  by  applying  separate  regression  coefficients  to 
events  below  and  above  the  water  table,  e.g., 

log  Yjbwl  =  a,bwt  +  b,bwt  (log  A0)jbwl  (5a) 

log  Y/WI  =  a r  +  b r  (log  A0)/Wt  (5b) 

Since  the  coupling  is  generally  significantly  smaller  for  events  above  the  water  table,  it  is 
likely  that  the  second  set  of  estimated  yields  (equations  5a  and  5b)  would  be  more  precise, 

since  different  relationships  for  events  below  and  above  the  water  table  are  used,  thus 

accounting  for  the  coupling  differences.  When  separate  coefficients  for  below  and  above  WT 
shots  are  used,  the  SD  of  log  Y,  S  is  given  by 


S2  =  d|  s,2  +  d2  s22 
d,  +  d2 

where  s1(  s2  are  the  SD  for  below  and  above  WT  shots  and  dj,  d2 
degrees  of  freedom  (see  Table  2). 


(6) 


are  the  corresponding 


Consider,  as  an  example,  plots  of  estimated  yield  based  on  A0.  When  the  regression 
coefficients  for  the  entire  data  set  are  used,  the  yield  estimates  for  events  below  the  water 
table  are  positively  biased  whereas  the  estimates  for  events  above  the  water  table  are  nega¬ 
tively  biased  and  the  SD  of  the  estimated  yield  is  0.276  (Table  2).  However,  when  separate 
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sets  of  coefficients  are  used,  the  systematic  biases  are  eliminated  and,  using  equation  (6),  the 


SD  reduces  to  0.196.  Furthermore,  the  correlation  coefficient  of  estimated  yield  with  known 
yield  improves  from  0.886  to  0.950  and  the  slope  increases  from  0.785  to  0.902.  A  similar 
improvement  is  obtained  for  yield  estimated  from  A/T,  A0  and  fc,  and  A/T  and  fc.  These 
improvements  are  summarized  in  Table  3  below. 


TABLE  3 

COMPARISON  OF  REGRESSION  RESULTS 


Slope 

Correlation  Coefficient 

SD  of 

log  Y 

Estimation 

Parameters 

Same 

Coefficients 

Separate 

Coefficients 

Same 

Coefficients 

Separate 

Coefficients 

Same 

Coefficients 

Separate 

Coefficients 

Ao 

0.785 

0.902 

0.886 

0.950 

0.276 

0.196 

A/T 

0.838 

0.941 

0.916 

0.970 

0.239 

0.152 

Ao.  L 

0.787 

0.964 

0.887 

0.982 

0.282 

0.126 

A/T,  fc 

0.866 

0.958 

0.930 

0.979 

0.224 

0.137 

At  least  some  of  this  increase  in  slope  and  correlation  coefficient  is  probably  due  to  the 
additional  degrees  of  freedom,  i.e. ,  the  use  of  twice  as  many  regression  coefficients  in  equa¬ 
tion  (5a)  or  (5b)  as  in  equation  (4).  Nevertheless,  a  comparison  of  the  SD  values  in  Table  3 
shows  that  substantially  better  yield  estimates  can  be  obtained  if  one  can  determine  whether 
the  event  occurred  above  or  below  the  water  table. 

Using  KEARSARGE  as  an  example,  the  derived  values  of  A0  and  fc  (indicated  on  Fig¬ 
ure  4)  and  the  regression  coefficients  in  Table  2  for  above  WT  shots  lead  to  the  estimated 
yield  value,  Y  given  by 


log  Y  =  -  0.211  +  0.915  (log  191.78)  +  0.374  (log  1.9) 


(7) 


which  gives  a  value  of  96  kt.  The  use  of  only  one  independent  variable  (viz.  A0)  will  provide 
a  value  of  77  kt.  Similarly,  use  of  regression  coefficients  for  all  data  lead  to  yield  estimates  of 
50  and  54  kt  when  results  with  one  (viz.  A0)  and  two  (viz.  A0  and  fc)  independent  variables, 
respectively  are  used.  As  expected,  the  first  value  (96  kt)  is  a  better  estimate  than  the  other 
three  values. 
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RESULTS  FROM  THE  RSTN  STATION,  RSSD 


Short-period,  vertical  component,  digital  data  from  17  NTS  explosions,  well  recorded  at 
the  RSTN  station,  RSSD  (40  samples/sec),  with  epicentral  distances  of  about  1300  km,  are 
also  analyzed  in  the  same  manner  as  the  JAS  data.  These  shots  are  listed  in  Table  4;  note  that 
10  shots  are  common  with  Table  1. 


TABLE  4 

NTS  SHOTS  RECORDED  AT  RSSD 


No. 

Date 

Name 

Lat 

Lon 

shot 

medium 

mb 

symbol 

I 

1982028 

28  jan  82 

JORNADA 

37.09 

-116.05 

tuff 

5.9 

X 

2 

1982175 

24  jun  82 

NEBBIOLO 

37.24 

-116.37 

rhyolite 

5.6 

o 

3 

1982217 

05  aug  82 

ATR1SCO 

37.08 

-116.01 

tuff 

5.7 

X 

4 

1982344 

10  dec  82 

MANTECA 

37.03 

-1 16.07 

alluvium 

4.6 

o 

5 

1983104 

14  apr  83 

TORQUOISE 

37.07 

-1 16.05 

tuff 

5.7 

X 

6 

1983125 

05  may  83 

CROWD1E 

37.01 

-1 16.09 

alluvium 

4.5 

0 

7 

1983146 

26  may  83 

FAHADA 

37.10 

-116.01 

tuff 

4.4 

0 

8 

1983160 

09  jun  83 

DANABLU 

37.16 

-116.09 

alluvium 

4.5 

0 

9 

1984091 

31  mar  84 

AGRIN1 

37.15 

-1 16.08 

alluvium 

4.1 

0 

10 

1984152 

31  may  84 

CAi  ruXTK 

3/.  10 

-1 16.05 

tuft 

5.8 

X 

1 1 

1984172 

20  jun  84 

DUORO 

37.00 

-116.04 

tuff 

4.6 

o 

12 

1984215 

02  aug  84 

CORREO 

37.02 

-116.01 

tuff 

4.7 

o 

13 

1985074 

15  mar  85 

VAUGHN 

37.06 

-116.05 

tuff 

4.8 

o 

14 

1985122 

02  may  85 

TOWANDA 

37.25 

-116.33 

tuff 

5.7 

X 

15 

1985163 

1 2  jun  85 

SALUT 

37.25 

-116.49 

rhyolite 

5.5 

X 

16 

1985339 

05  dec  85 

KINIBOTO 

37.05 

-116.05 

tuff 

5.7 

X 

17 

1986318 

14  nov  86 

GASCON 

37.10 

-116.05 

tuff 

5.8 

X 

* 

O  represents 

above  and  X 

below  the  water 

table 

Spectra  of  Bn  (window  length  6.4  sec)  from  ATR1SCO  and  NLBBIOLO.  two  shots  of 
similar  magnitudes  but  the  first  below  WT  and  the  other  above  WT,  are  shown  in  Figure  S. 
The  spectra  are  not  smoothed  and  arc  not  corrected  for  instrumental  response.  The  S/N  is 
gv  >d  up  to  about  5  Hz.  The  attenuation  correction  for  the  RSSD  data  is  assumed  to  be  t*  ■- 
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Figure  8.  Venical-component  displacement  amplitude  spectra  of  Pn  (6.4  sec  long  window), 
similar  to  those  in  Figure  2,  for  ATRISCO  (below  water  table)  and  NEBBIOLO  (above  water 
table)  recorded  at  RSSD.  S/N  appears  good  up  to  about  5  Hz  only. 


0.4  as  the  source-receiver  distance  is  about  4  times  that  for  the  station  JAS.  The  correspond¬ 
ing  displacement  spectral  fits  and  the  derived  values  of  fc  and  A0  are  shown  in  Figure  9.  Note 
that  fc  for  the  above  WT  shot  is  considerably  higher  than  for  the  below  YVT  shot,  similar  to 
the  result  in  Figure  4. 

A  plot  of  ft  versus  A0  for  all  17  shots  recorded  at  RSSD  is  shown  in  Figure  10a;  the 
explosion  numbers  correspond  to  those  in  Table  4.  The  mean  slope  and  intercept  values  are 
-0.398  (0.038)  and  1.071  (0.090),  respectively.  When  the  mean  slope  derived  from  the  entire 
data  set  is  forced  on  the  below  WT  data  alone  (Figure  10b),  the  mean  intercept  value  is  1.065 
(0.032)  whereas  the  corresponding  result  for  only  above  WT  shots  (Figure  10c)  is  1.076 
(0.044).  Thus  there  is  hardly  any  separation  between  the  below  and  above  WT  populations. 

Time-domain  measurements  of  log  A/T,  a  measure  of  inb,  are  made  in  the  same  manner 
as  for  the  JAS  recordings.  Plots  of  fc  versus  log  A/T  are  shown  in  Figures  11a.  b,  and  c  for 
all  (17)  shots,  8  below  WT  shots,  and  9  above  WT  shots,  respectively.  Forcing  the  mean  slope 
of  -0.480  obtained  from  the  entire  data  set  (Figure  11a),  the  intercept  values  for  below  and 
above  WT  populations  are  1.102  (0.058)  and  1.189  (0.057),  respectively.  Therefore,  the  comer 
frequencies  for  a  given  mb  for  shots  below  and  above  WT  are  not  significantly  d  fferent. 

A  plot  of  A()  versus  yield  (on  log-log  scale)  for  all  17  shots  indicated  the  mean  slope  of 
1.648,  considerably  larger  than  the  theoretically  expected  value  of  1  or  the  observed  value  of 
close  to  1  for  the  JAS  data.  One  reason  for  this  is  the  lack  of  sufficient  overlap  in  the  yield 
values  of  below  and  above  WT  shots.  Forcing  a  slope  of  1  to  the  data  punts,  the  mean  inter¬ 
cept  values  for  all  (17)  shots,  X  below  WT,  and  9  above  WT  shots  were  0.596  (0  103),  0.954 
(0  1  14),  and  0.278  (0  OsXi,  respectively.  The  below  and  above  WT  populations  are  therefore 
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Figure  10a.  fc  versus  A„  for  8  below  water  table  (denoted  by  X)  and  9  above  water  table 
(marked  by  circles)  shots  numbered  as  in  Table  4  and  based  on  t*  -  0.4.  The  least  squares 
linear  regression  through  all  17  data  points  (continuous  line)  has  correlation  coefficient  of 
0.930,  slope  of  -0.398,  and  intercept  of  1.071.  One  SI)  variation  in  the  intercept  value  is 
shown  by  the  dashed  lines 
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rigure  1  la.  fc  versus  log  A/T  for  8  below  water  table  (denoted  by  X)  and  9  above  water  table 
(marked  by  circles)  shots  based  on  t*  =  0.4.  The  least  squares  linear  regression  through  all  17 
data  points  (continuous  line)  has  correlation  coefficient  of  0.851,  slope  of  0  480,  and  intercept 
of  1. 148.  One  SU  variation  in  the  intercept  value  is  shown  bv  the  dashed  lines. 


separated  by  0.676  (0.128).  This  is  close  to  the  separation  observed  from  22  shots  recorded  at 
JAS. 
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COMPARISON  WITH  RESULTS  FROM  LARGER  DATA  BASES 


Blandford  and  Klouda  (1980)  measured  regional  phases  including  Pn  from  a  large 
number  of  NTS  explosions  recorded  at  TFO,  at  epicentral  distances  of  about  550  km.  Their 
results  for  magnitudes  based  on  amplitude  of  Pn  (the  "b"  phase)  are  used  here  to  compare  the 
magnitude  versus  yield  relationships  for  below  above  WT  shots.  The  Rainier  Mesa  shots  are 
considered  to  be  below  WT  because  of  a  perched  water  table  at  shallow  depths  in  the  region. 
A  plot  of  magnitude  versus  log  yield  for  all  available  below  WT  (35)  shots  shows  a  correla¬ 
tion  coefficient  of  0.918,  mean  slope  of  0.744  (0.056)  and  mean  intercept  of  4.477  (0.104). 
Corresponding  results  from  above  WT  (112)  shots  are:  correlation  coefficient  0.796,  mean 
slope  0.899  (0.065),  and  mean  intercept  3.685  (0.065).  These  plots  indicate  a  magnitude 
difference  of  about  0.46  between  below  and  above  WT  shots  of  150  kt  each.  Note  that  this 
value  is  fairly  close  to  the  separation  of  0.36  (0.05)  obtained  for  shots  recorded  at  JAS. 

Using  amplitude  data  from  the  ISC  bulletins,  P.  D.  Marshall  (written  communication) 
obtained  maximum  likelihood  estimates  of  the  seismic  magnitude  mb  of  a  large  number  of 
NTS  explosions.  The  methodology  used  is  described  in  detail  in  Marshall  ct  al.  (1986).  Out 
of  a  total  of  300  events,  287  events  were  found  suitable  for  studying  the  effect  of  the  water 
table.  Plots  of  mb  versus  log  yield  for  below  WT  (114)  shots  indicate  correlation  coefficient  of 
0.971,  mean  slope  of  0.754  (0.018),  and  mean  intercept  of  4.096  (0.035).  Results  from  173 
above  WT  shots  are:  correlation  coefficient  0.872,  mean  slope  0.899  (0.039),  and  mean  inter¬ 
cept  3.491  (0.054).  These  plots  will  indicate  a  magnitude  difference  of  about  0.29  between 
below  and  above  WT  shots  of  150  kt  each. 
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An  attempt  is  now  made  to  estimate  the  effect  of  water  table  when,  considering  the  size 
of  the  explosion  cavity  and  the  surrounding  fracture  zone,  shots  may  be  considered  to  be  well 
separated  from  the  water  table.  A  rough  measure  of  the  explosion  cavity  radius  Rc  is  given 
by 

Rc=10Yin  (8) 

where  Y  is  the  yield  in  kt.  This  relationship  is  not  much  different  from  the  empirical  depth- 
and  medium-dependent  relationships  for  NTS  explosion  cavity  radii  obtained  by  Closmann 
(1969).  One  may  consider  shots  having  their  shot  points  at  least  2  Rc  below  or  above  WT  to 
be  well  separated  from  the  water  table.  mb  versus  log  yield  plot  for  72  well  below  WT  shots 
(including  Rainier  Mesa  shots  which  are  all  assumed  to  be  well  below  WT)  indicated  correla¬ 
tion  coefficient  of  0.981,  mean  slope  of  0.746  (0.017),  and  mean  intercept  of  4.118  (0.035). 
The  correlation  coefficient  for  the  data  from  125  well  above  WT  shots  is  0.708,  mean  slope  is 
0.711  (0.064),  and  mean  intercept  is  3.670  (0.075).  These  plots  indicate  the  magnitude  of  a 
below  WT  shot  to  be  larger  than  that  of  an  above  WT  shot  by  about  0.52,  if  each  has  the 
same  yield  of  150  kt.  As  expected,  the  separation  between  the  below  and  above  populations 
is  enhanced  when  only  those  shots  that  are  well  below  and  well  above  WT  are  considered. 
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EFFECT  OF  C.AS  POROSITY  AND  CAVITY  RADII  FOR  SHOTS  IN  TUFF 

Magnitude  versus  yield  plots  for  above  WT  shots  show  considerably  more  scatter  than 
those  for  below  WT  shots.  According  to  Patton  (1988),  a  possible  reason  for  this  is  the  effect 
of  gas-filled  voids  in  the  emplacement  medium  on  source  coupling.  His  investip"‘;o"  b°sed  on 
the  use  of  mb(Lg)  magnitudes  showed  the  gas  porosity  to  significantly  reduce  mb  (Lg)  for 
shots  in  tuff  (and  rhyolite  since  no  distinction  was  made  between  the  two  media,  H.  J.  Patton, 
oral  communication).  He  also  examined  the  effect  of  gas  porosity  on  Marshall’s  mb  values 
for  21  Yucca  shots  and  found  the  effect  to  be  very  similar  to  that  for  mb  (Lg).  In  this  study, 
we  examine  a  much  larger  data  set  consisting  of  a  total  of  148  shots  (91  below  and  57  above 
WT)  from  the  Yucca  Flats,  Pahute  Mesa,  and  Rainier  Mesa  regions  of  the  NTS.  These  are 
the  only  tuff-rhyolite  shots  for  which  gas  porosity  values  are  available,  out  of  a  total  of  287 
shots  for  which  Marshall’s  mb  are  known. 

Approximate  expressions  for  cavity  radii  in  tuff  and  rhyolite,  based  on  Closmann’s 
(1969)  study  (see  also  Mueller  and  Murphy,  1971),  are 

Rc  (tuff)  =  31.41  y029  h-011  (9) 

and 

Rc  (rhyolite)  =  28.66  Y029  h'°  1 1  (10) 

where  h  is  the  shot  depth  in  meters.  These  two  equations  were  used  to  determine  cavity  radii 
for  all  148  tuff-rhyolite  shots.  Explosions  with  their  shot  points  at  least  two  cavity  radii 
below  the  water  table  and  all  Rainier  Mesa  shots  may  be  considered  to  be  well  below  ihc 
water  table.  There  were  56  such  shots  with  known  gas  porosity  and  8  more  for  which  porosity 
was  not  known.  For  these  56  +  8  =  64  well  below  the  water  table  shots,  mb  versus  yield  rela 
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tionship  was  found  to  be 


mb  =  0.744  (0.018)  logY  +  4.117  (0.036).  (11) 

The  SD  of  estimated  mb  is  0.109. 

A  jdA  of  mb  residuals  (obsep/ed  -  f’S!:rr,ated\  where  the*  mated  values  are  derived  for 
known  Y  from  equation  (11),  versus  gas  porosity  (in  volume  %)  is  shown  in  Figure  12a  for 
all  148  shots.  The  usual  least  squares  regression  (continuous  line),  based  on  die  assumption 
that  porosity  is  error-free,  has  a  mean  slope  of  -0.033  (0.002)  and  the  SD  of  mb  residual  is 
0.13.  Note  that  several  porosity  values  in  Figure  12a  are  negative  "even  though  physically 
unreal  as  a  result  of  small  errors  in  the  measured  parameters  input  to  the  calculation" 
(Howard,  1983).  Furthermore,  there  are  a  large  number  of  shots  with  exactly  zero  porosity. 
All  this  means  that  the  error  in  porosity  is  substantial  and  should  not  be  neglected.  In  order 
to  study  the  effect  of  errors  in  both  x  and  y  variables,  where  x  denotes  porosity  and  y 
represents  mb  residual,  the  maximum  likelihood  linear  fitting  method  of  Ericsson  (1971)  is 
used  to  obtain  mean  slope  and  intercept  values.  The  variables  x  and  y  are  assumed  to  be  com¬ 
pletely  uncorrelated  and  the  ratio  g  h  oy/ox,  where  o  denotes  the  SD,  is  assumed  to  be 
known  Since  the  smallest  value  of  observed  porosity  for  several  shots  is  -5,  a  rough  estimate 
of  ax  may  be  assumed  to  be  5.  Similarly,  from  the  data  used  in  equation  (11).  av  is  about  0.1 
so  that  g  may  be  taken  to  be  0.02.  The  resulting  maximum  likelihood  fit,  with  derived  mean 
slope  value  of  -0.042  (0.003),  is  also  shown  in  Figure  1 2a  (dashed  line).  When  56  well  below 
WT  shots  are  excluded,  the  results  for  the  remaining  92  shots  are  shown  in  Figure  12b.  When 
porosity  (x)  is  assumed  to  be  error-free;  the  mean  slope  is  -0.034  (0.002)  whereas  the  value 
increases  to  -0.043  (0.003)  if  g  0.02.  If  only  shots  well  above  WT  are  considered,  the 
regressions  for  31  such  shots  is  shown  in  Figure  12c.  Hie  mean  slope  changes  from  -0.026 


41 


MB  [MARSHALL]  RESIDUAL 

.20  -0.80  -0.M0  -0.03  0.  HO 


'-5.00  0.00  5.00  10.00  15.00  20.00  25.00 

GAS  POROSITY  (/  VOLUME] 

Figure  12a.  mb  residual  versus  gas  porosity  (in  volume  %)  for  148  tuff-rhyolite  shots  with 
known  gas  porosity  values.  Below  and  above  water  table  shots  are  denoted  by  X  and  circle, 
respectively.  The  least  squares  linear  regression  (continuous  line)  has  a  mean  slope  of  -0.033 
(0.002)  whereas  regression  based  on  g  =  0.02  (dashed  line)  shows  a  larger  slope  of  -0.042 
(0.003). 
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Figure  12b.  Similar  to  Figure  12a  for  92  shots  obtained  by  excluding  56  well  below  water 
table  shots  from  those  used  in  Figure  12a.  The  least  squares  linear  regression  (continuous 
line)  has  a  mean  slope  of  -0.034  (0.002)  whereas  regression  based  on  g  =  0.02  (dashed  line) 
shows  a  larger  slope  of  -0.043  (0.003). 
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Figure  12c.  Similar  to  Figure  12a  for  31  well  above  water  tab'e  shots.  The  least  squares 
linear  regression  (continuous  line)  has  a  mean  slope  of  -0.026  (0.006)  vhereas  regression 
based  on  g  =  0.02  (dashed  line)  shows  a  larger  slope  of  -0.054  (0.016). 
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(0.006)  for  error-free  x  to  -0.054  (0.016)  for  g  =  0.02.  Note  that  the  mean  slopes  values  based 
on  porosity  not  being  error-free  (dashed  lines)  are  all  significantly  larger  than  those  based  on 
error-free  x  (continuous  lines),  thereby  suggesting  enhanced  dependence  of  mb  residual  on  gas 
porosity.  All  three  figures  indicate  gas  porosity  to  reduce  the  observed  m5  by  significant 
amounts. 

We  now  consider  the  effect  of  gas  porosity  on  shots  close  to  the  water  table,  i.e.  on 
92-31  =61  shots  that  are  neither  well  below  nor  well  above  WT  (i.e.  shots  within  2  Rc  of 
WT).  Plots  of  mb  residual  versus  porosity  from  all  61  shots,  a  subset  of  35  that  are  below 
WT,  and  the  remaining  26  above  WT  shots  are  shown  in  Figures  13a,  13b,  and  13c,  respec¬ 
tively.  The  mean  slopes  are  -0.029  (0.003),  -0.013  (0.007),  and  -0.027  (0.004)  based  on 
error-free  x,  and  -0.039  (0.004),  -0.105  (0.025),  and  -0.036  (0.007)  for  g  -  0.02,  respectively. 
All  slope  values  are  negative,  confirming  the  expected  effect  of  gas  porosity.  Since  these  61 
shots  are  all  near  WT,  they  should  also  be  expected  to  show  the  effect  of  proximity  of  their 
cavity  radii  and  fracture  zone  to  the  water  table.  Defining  "scaled  distance  to  water  table" 
(SDWT)  as  the  vertical  distance  to  the  water  table  (considered  positive  for  below  and  negative 
for  above  WT  shots,  respectively)  from  the  shot  point  divided  by  Rc  ,  where  Rc  is  given  by 
equations  (9)  or  (10),  the  results  for  61  (35  below'  and  26  above  WT)  shots  are  shown  in  Fig¬ 
ure  14.  Linear  regression  of  the  data  (continuous  line)  has  a  mean  slope  of  0.141  (0.017),  sug¬ 
gesting  that  explosion  cavity  and  its  proximity  to  the  water  table  significantly  influence  the 
source  coupling.  The  lack  of  a  sharp  discontinuity  where  the  two  populations  meet  (i.e.  at 
SDWT  -  0)  means  that  the  distinction  between  below  and  above  WT  shots  is  lost  for  shots 
very  close  to  WT.  For  explosions  with  their  shot  points  below  WT,  shot-generated  cavity  and 
surrounding  fractures  extending  above  the  water  table  (especially  for  shots  with  their  cavities 
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Figure  13a.  mb  residual  versus  gas  porosity  (in  volume  %)  for  61  tuff-rhyolite  shots  within  2 
cavity  radii  ot  the  water  table.  Below  and  above  water  table  shots  are  denoted  by  X  and  cir¬ 
cle,  respectively.  The  least  squares  linear  regression  (continuous  line)  has  a  mean  slope  of 
-0.029  (0.003)  whereas  regression  based  on  g  =  0.02  (dashed  line)  shows  a  larger  slope  of 
-0.039  (0.004). 
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Figure  13b.  Similar  to  Figure  13a  for  35  tuff-rhyolite  shots  that  are  below  but  within  two  eas¬ 
ily  radii  of  the  water  table.  The  least  squares  linear  regression  (continuous  line)  h  is  a  mean 
slope  of  -0.013  (0.007)  whereas  regression  based  on  g  =  0.02  (dashed  line)  shows  a  larger 
slope  of  -0.105  (0.025). 
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Figure  13c.  Similar  to  Figure  13a  for  26  tuff-rhyolite  shots  that  are  above  but  within  two  cav¬ 
ity  radii  of  the  water  table.  The  least  squares  linear  regression  (continuous  line)  has  a  mean 
slope  of  -0.027  (0.004)  whereas  regression  based  on  g  =  0.02  (dashed  line)  shows  a  larger 
slope  of  -0.036  tO.007). 
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Figure  14.  mb  residual  versus  scaled  distance  to  water  table  (SDWT)  for  61  tuff-rhyolite  shots 
(35  below  and  26  above  WT,  denoted  by  X  and  circle,  respectively)  within  two  cavity  radii  of 
W  F.  No  sharp  discontinuity  is  observed  between  the  below  and  above  WT  populations  (at 
SDWT  =  0).  The  least  squares  regression  (continuous  line)  has  a  mean  slope  of  0.141  (0.017) 
and  the  SD  of  inh  residual  is  0.135. 


crossing  the  water  table)  appear  to  reduce  the  coupling.  Similarly,  for  above  WT  shots,  explo¬ 
sion  cavity  and  fractures  intersecting  the  water  table  seem  to  enhance  the  coupling. 
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DISCRIMINATION  BETWEEN  BELOW  AND  ABOVE  WT  SHOTS 


Analysis  of  the  Pn  spectra  of  22  shots  recorded  at  JAS  shows  some  separation  on  plots 
of  fc  versus  A0  derived  from  the  Pn  spectra  of  22  shots  recorded  at  JAS.  However,  the 
separation  is  probably  not  large  and  significant  enough  to  be  of  much  practical  value  (see  Fig¬ 
ures  5  and  6).  On  the  other  hand,  plots  of  A0  versus  yield  show  the  largest  separation  between 
the  two  types  of  shots.  It  seems  therefore  that  the  low- frequency  spectral  amplitude  of  Pn  is  a 
possible  discriminant.  This  is  also  supported  by  the  spectral  ratio  of  KKARSARGE  (above 
WT  shot)  to  JORNADA  (below  WT  shot)  in  Figure  3  in  which,  at  low  frequencies,  the  ampli¬ 
tude  ratio  is  much  smaller  than  1. 

The  low-frequency  amplitude  values  may  also  be  obtained  by  bandpassing  the  data  and 
computing  the  root-mean-square  (RMS)  amplitudes,  after  correcting  for  the  instrumental 
response  and  noise  level  before  the  onset  of  direct  Pn.  Results  for  the  05-1.5  Hz  passband, 
with  window  length  of  6.4  see,  for  the  22  shots  recorded  at  JAS  are  shown  in  Figuie  i5.  the 
least  squares  regression  lines  for  the  two  types  of  shots,  assuming  a  slope  of  1,  are  separated 
by  0.542  magnitude  unit.  These  results  are,  as  expected,  nearly  identical  with  those  tor  A,, 
The  below  WT  shots  show  larger  deviations  from  the  mean  than  the  above  W  I  shots.  The 
RMS  value  from  one  shot  (SAL(JT)  in  Figure  15.  denoted  by  the  letter  P,  is  close  to  that 
expected  for  an  above  WT  shot.  A  possible  reason  for  this  may  be  the  fact  that  this  is  the 
only  below  WT  shot  with  non-7.ero  gas  porosity  value  (viz.  4%). 

An  impoitant  factor  contributing  to  the  scatter  in  Figure  15  is  perhaps  proximity  <>!  the 
shot  cavities  of  several  explosions  to  the  water  table  Out  of  a  total  of  22  shots,  the  shot 
p. . mts  of  3  are  within  one  cavity  radius  of  W|  (marked  by  long  arrows)  and  those  of  d  other 


LOG  YIELD 

f  igure  15.  RMS  amplitudes  of  Rn.  for  the  bandpass  0  5  to  1.5  Hz,  versus  yield  for  22  shots 
recorded  at  JAS.  The  mean  lines  with  slope  of  unity  for  the  below  and  above  WT  populations 
are  separated  bv  0.542  magnitude  unit  The  arrows  indicate  explosions  with  shot  points  within 
two  cavitv  radii  of  the  water  table. 


are  within  two  cavity  radii  of  WT  (denoted  by  small  arrows).  If  this  were  not  the  case,  the 
RMS  values  for  these  12  shots  would  move  towards  the  direction  pointed  by  the  arrows, 
thereby  resulting  in  somewhat  larger  separation  between  the  two  types  of  shots. 

A  plot  of  the  RMS  values  for  Lg  with  window  length  of  51.2  sec  (but  otherwise  obtained 
in  the  same  manner  as  for  Pn)  versus  yield  for  10  shots  for  which  Lg  was  available  (most  oth¬ 
ers  were  clipped)  is  showm  in  Figure  16.  The  below  and  above  WT  populations  are  separated 
by  0.300  magnitude  unit.  It  seems  therefore  that  Lg  is  considerably  less  sensitive  than  Pn  to 
whether  the  shot  is  below  or  above  WT  and  Pn/Lg  amplitude  ratio  is  a  possible  discriminant. 
This  is  in  agreement  with  an  analysis  of  Blandford  and  Klouda's  (1980)  data  showing  larger 
separation  between  below  and  above  WT  shots  for  mh  (Pn)  than  for  mb  (Lg). 


DISCUSSION  AND  CONCLUSION 


Records  of  JORNADA  and  KEARSARGE,  with  good  S/N,  are  also  available  at  the  SRO 
station,  ANMO  at  epicentral  distance  of  about  900  km.  The  two  records  are  shown  in  Figure 
17.  Note  that  just  as  in  Figure  1,  JORNADA  has  much  larger  amplitudes  than  KEARSARGE 
in  the  first  several  seconds  of  Pn  but  not  so  in  the  later  wavetrain.  The  spectra  of  Pn  (6.4  see 
window)  for  the  two  shots  are  shown  in  Figure  18;  S/N  is  good  up  to  about  6  Hz.  The  spec¬ 
tral  ratio  KEARSARGE/JORN ADA  is  shown  in  Figure  19  and  appears  similar  to  that  in  Fig¬ 
ure  3  although  the  mean  slope  value  is  slightly  larger.  The  spectral  fits,  based  on  t*  =  0.25, 
are  shown  in  Figure  20  and  show  differences  between  the  two  shots  that  are  similar  to  those 
in  Figure  4.  One  may  thus  conclude  that  the  differences  in  the  spectral  characteristics  of  these 
two  shots,  as  observed  at  ANMO,  are  similar  to  those  observed  at  JAS  and  the  nearby  station 
C.MB. 

A  comparison  of  the  spectral  characteristics  of  below  and  above  WT  shots  shows  that  for 
most  parameters,  there  is  a  lack  of  systematic  differences  and  considerable  overlap  between 
the  two  types  of  shots,  even  though  records  from  the  same  station  are  used  so  that  source- 
receiver  paths  are  not  much  different.  This  is  not  surprising  in  view  of  the  extreme  variability 
of  regional  phases  from  NTS  explosions  at  essentially  the  same  location  noted  by  Springer 
and  Denny  (1976,  especially  their  Figure  4)  and  Gupta  cr  al  (  1984,  see  their  Figure  1  1).  Plots 
of  fc  versus  A0  appear  to  show  some  separation  between  below  and  above  WT  shots  recorded 
at  JAS  (Figures  5b.  c  and  6b,  c).  This  separation  is,  however,  not  evident  in  the  RSSD  data 
(Figures  10b,  c),  perhaps  because  the  RSSD  results  are  not  as  reliable  due  to  poor  overlap  in 
the  fc  and  A(j  values  of  the  below  and  above  WT  populations.  Moie  data  from  several  other 
stations  should  be  examined  to  explore  further  the  possibility  of  distinguishing  between  below 
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Figure  17.  ~  -period,  vertical -component  records  of  JORNADA  (below  water  table)  and 
KEARSARGE  (above  water  table)  at  ANMO,  similar  to  those  in  Figure  1.  Note  that  in  the 
first  several  seconds  of  Pn,  KEARSARGE  has  much  smaller  amplitude  than  JORNADA. 
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Figure  19.  Spectral  ratio  KEARSARGE/JORNADA  of  Pn  at  ANMO,  similar  to  that  in 
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Figure  3. 
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and  above  WT  shots  from  measurements  of  their  fc  and  A0  values.  Broadband  instruments  are 
recommended  so  that  both  A0  and  fc  are  better  constrained.  A  determination  from  seismic  data 
alone  of  whether  a  shot  was  fired  below  or  above  WT  is  important  since  it  can  substantially 
improve  yield  estimates.  Plots  of  ft  versus  A/T  and  fc  versus  yield  show  no  observable  dis¬ 
tinction  between  the  two  types  of  shots.  The  largest  systematic  difference  between  below  and 
above  WT  shots  lies  in  plots  of  A0  versus  yield  and,  to  a  somewhat  smaller  degree,  in  plots  of 
A/T  (a  measure  of  mb)  versus  yield. 

Analysis  of  both  regional  and  teleseismic  mb  data  from  substantially  larger  data  bases 
also  demonstrate  that,  for  a  given  yield,  mb  values  for  shots  below  WT  are  generally  larger 
than  those  for  shots  above  WT  by  0.3  to  0.5  magnitude  unit.  Gas  porosity  of  the  shot  medium 
appears  to  play  an  important  role  in  the  determination  of  this  large  magnitude  bias.  Yield  esti¬ 
mates  can  be  significantly  improved  if  it  is  known  whether  the  shot  point  is  below  or  above 
the  wu'er  table.  A  knowledge  of  the  gas-filled  porosity  of  the  shot  medium  may  further 
improve  the  precision  of  yield  determinations  of  shots  above  the  water  table. 

A  surprising  and  interesting  result  is  the  important  role  played  by  explosion  cavity  radii 
in  determining  the  source  coupling.  The  correlation  between  mb  residual  and  SDWT  (positive 
for  below  and  negative  for  above  WT  shots)  for  shot  points  within  two  cavity  radii  of  WT 
makes  a  strong  case  for  the  intluence  ot  explosion  cavity  and  associated  fracture  zone  on 
source  coupling.  Significant  improvement  in  yield  estimates  can  be  achieved  by  combining 
mb  values  with  both  gas  porosity  and  SDWT  where  the  latter  may  be  computed  using  empiri¬ 
cal  medium  and  depth-dependent  relationships  lor  explosion  cavity  radii.  A  comparison  ot 
1’n  and  l.e  pluses  recorded  at  a  common  station  may  be  eilective  in  distinguishing  below  and 
above  WT  shots  except  those  that  are  fired  close  to  the  water  table. 
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